set more off


local J `1'

forvalues i = 1/`J' {

gen epsilon_`i' = -ln(-ln(runiform())) 
}
